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Abstract 

A method is given for creating material with a desired refraction co- 
efficient. The method consists of embedding into a material with known 
refraction coefficient many small particles of size a. The number of parti- 
cles per unit volume around any point is prescribed, the distance between 
neighboring particles is 0(a a ) as a ^ 0, < k < 1 is a fixed parame- 
ter. The total number of the embedded particle is 0(a K ~ 2 ). The physical 
properties of the particles are described by the boundary impedance £ m 
of the m — th particle, £ m = 0(a~ K ) as a — » 0. The refraction coefficient 
is the coefficient n 2 {x) in the wave equation [V 2 + k 2 n 2 (x)]u — 0. 



PACS: 03.4Q.Kf, 03.50 De, 41.20.Jb, 71.36. +c 
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1 Introduction 

The problem we are concerned with is the following: how does one create in a 
given bounded domain D c R 3 a material with a desired refraction coefficient 
n 2 (x). The domain D originally is assumed to be filled with a material with a 
known refraction coefficient n^(x). We assume that Im n^(x) > and n^(x) = 1 
in D' := R 3 \ D. Originally the wave equation is: 

L Q u := [V 2 + k 2 nl(x)]u = in R 3 , k = const > 0, (1) 

u = e ika - x + v , aeS 2 , (2) 
e lkr /1\ , x 



v =A (p,a,k) ho - , r = |ar|->oo, /3 := -. (3) 

r \r J r 

The function vq is the scattered field, Ao(/3,a,k) is the scattering amplitude, 
u(x, a, k) is the scattering solution, S 2 is the unit sphere in R 3 . 
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We embed M small particles D m , S m := dD m , 1 < m < M, into D, so that in 
any subdomain Acfl there are 

Af(A) = -J— / N(x)dx[l +o(l)], a^O (4) 



small particles. Here N(x) > is a continuous (Or piecewise-continuous) func- 
tion which we can choose as we wish, < k < 1 is a parameter which is our 
disposal. For simplicity we assume that particles D m are balls centered at the 
points x m and of radius a independent of m. The distance d between neighbor- 
ing particles is assumed to be 

d = 0{a^) asa^O. (5) 

The properties of a particle are described by the boundary impedance 

Cm = (6) 

where h(x) is a continuous function on D, Im h(x) < 0. The function h(x), 
as N(x), we can choose as we wish. The scattering solution u(x,a,k) in the 
presence of the embedded particles solves the problem: 

L u = in M 3 \ U^f =1 D m , (7) 

u N + CmU = on S m , 1 < m < M, (8) 
u = uo(x, a, k) + v, (9) 

gifcr / t \ x 

v = Ai((3,a,k) ho - , |ar|=r-^oo, /?:=-. (10) 

r \r J r 

Let us now describe our results. We prove that problem (fT|)- (fiT))) has a unique 
solution u(x, a, k) :— Uja(x,a,k). We prove that given an arbitrary function 
n 2 (x) such that n 2 (x) = 1 in D' , n 2 (x) is continuous or piecewise-continuous 
in D (with the set of discontinuities of Lebesgue measure zero in R 3 ), one can 
choose N(x) and h(x) so that the limit 

ip := ip(x, a, k) — lim %(i,a,fe) (11) 

M— *oo 

exists and satishes the equation 

[V 2 + k 2 n 2 {x)]%l) = inM 3 , (12) 
xp = uq{x, a, k) + w(x, a, k), (13) 

ip = e ikax +A(/3,a,k)—+o(-\ r = \x\ -> 00, (3 := -. (14) 
r \r J r 

Therefore the medium with embedded particles in the limit M — > 00, or, which 
is the same by ((4]), in the limit a — ► 0, has a desired refraction coefficient n 2 (x). 
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In Section 2 we formulate the recipe for choosing N(x) and h(x) which guaran- 
tees the existence of the limit (fTTj) which solves problem (fT2|) - (fT4|) . We do not 
assume that the small particles are embedded periodically. 
The aim of this paper is to make clear for a wide audience of engineers and 
physicists our recipe for creating material with any desired refraction coefficient 
and to formulate two technological problems which must be solved in order 
that our theory can be immediately implemented experimentally. Theoretical 
justification of our results are given in [T]-[3], see also [I]- [5J. 

2 The recipe for creating material with a desired 
refraction coefficient 

The problem we are interested in is the following: 

One is given tIq{x) and wants to create a refraction coefficient n 2 (x). 

Here is our recipe for doing this. 
Step 1. Calculate the function 

p(x) := k 2 [nl(x) ~ n 2 (x)} := pi(x) + ip 2 (x), (15) 
where pi(x) = Re p(x), p2 = Im p(x). 

Step 2. Find two functions N(x) > and h(x) = hi(x)+ili2(x) from the relation 

<brh(x)N(x) = p{x). (16) 

This can be done by infinitely many ways. For example, one may fix N(x) > 
and define 

hl = P1{X )„ h 2 = P2 ^\. (17) 

If one wishes to deal only with passive materials, then one requires Im n 2 (x) > 0, 
Im h(x) < 0, and, if Im it,q(x) < Im n 2 (x),then Im p(x) < 0. 
Step 3. Partition the domain D into a union of small cubes A p , 1 < p < P , 
without common interior points, D = Up =1 A p , the center of A p is denoted by 

y p , the side of A p is of the order O(a^). In each cube A p embed AA(A p ) small 
particles, where Af(A p ) is defined in ((4]). The distance d between neighboring 
particles should be d = 0(a~5~). The given order of the smallness of d as a — > 
is important, but the distance need not be exactly the same. The boundary 
impedance of each of the small particle embedded in A p make equal to h ^ y ^ ; 
where h(x) = hi(x) + ih,2(x) is the function found in Step 2 of the recipe. 

Theorem 2.1. After the completion of Step 3, the material, obtained from 
the original one with the refraction coefficient Uq(x), will have the refraction 
coefficient n 2 M {x), and YanM^oo n \i( x ) — n 2 (x). 

Proof of this theorem one finds in papers pQ and [3] . 
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3 A discussion of the recipe 



Step 1 of the recipe is trivial. Step 2 is also trivial. One may choose N(x) > 
to satisfy some practical requirements. For example, if one chooses N(x) small, 
then the total number of particles will be smaller. Practically one cannot take 
the limit M — > oo, i.e., in the limit a — ► 0, and one stops at some finite value 
of M, or of a > 0. The two technological problems, that have to be solved in 
order that our recipe can be implemented experimentally, are: 

1) How does one embed a small particle at a given point into the given material 
inD? 

2) How does one prepare a small particle, a ball of radius a centered at a point 
x m , with the prescribed boundary impedance £ m = MggJ ? 

Here h(x) is the function, found at Step 2 of the recipe? 

Possibly, the first technological problem can be solved by the stereolitography 
process. 

One should be able to solve the second technological problem because its 
limiting cases £ = (acoustically hard particles, particles from insulating mate- 
rial) and ( m — oo (acoustically soft particles, perfectly conducting particles) are 
easy to solve in practice, so the intermediate values of the boundary impedance 
should be also possible to prepare. 

A similar theory has been developed in paper [5] for electromagnetic wave scat- 
tering by many small dielectric and conducting particles embedded in an inho- 
mogeneous medium. 

4 Electromagnetic waves 

Assume now that the governing equations are the Maxwell equations 

VxE = iufiH, V x H = -iwe'{x)E in R 3 , (18) 

u = const, e'(x) = e = const in D 1 , lo > is frequency, e'(x) = e(x) + 
^77^7 o~{x) > is the conductivity, a(x) = in D' . We assume that e'(x) € 
C 2 (1R 3 ), e'(x) 7^ 0, is a twice continuously differentiable function. Let k = — , 
c = ui-y/ejl is the wave velocity in D'. The incident plane wave is £e lka x , a £ 
S 2 , a ■ £ = 0, £ is a constant vector. 

Under the above assumptions the electrical field E(x) is the unique solution 
to the equation (see [2]): 



Eo(x) = Se* ka -" + / g(x,y)p(y)E (y)dy 

+ V X / g{x,y)q(y) ■ E (y)dy, g{x,y) := 



(19) 



A-k\x — y\ ' 



where 



p(x):=K 2 (x)-k 2 , K\x):=u, 2 e'(x)n; q {x):=^^±. (20) 



4 



If M small particles D m , 1 < m < M, are embedded in D, then the basic 
equation (fl"9]l becomes 



M „ M 



E M {x) = E Q (x)+ V] / g(x,y)p{y)E M (y)dy+Y] / g{x,y)q{y)E M {y)dy. 

(21) 

It is proved in [5] that if the size a of small particles tends to zero, if the number 
of these particles in any open subset A of D is 

M(A) = -^ [ N(x)dx[l+o(l)}, a^O, (22) 

a J A 

and if the distance d between neighboring particles is d = 0(a S ~3~), then there 
exists the limit limM^o Em (x) = E e (x) . 

The limiting field E e (x), i.e., the effective field in the medium, solves the 
equation: 

E e (x)=E (x)+ ( g(x,y)C(y)E e (y)dy, C(y) = N(y)c(y), (23) 
Jd 

where 

c(y) = lim a 3 -" / p{x)dx. (24) 

J\y~x\<a 

If, e.g., the small particle D m is a ball of radius a centered at a point y, and 

2lvl{l-\A)\ \x\<a: 



P(X) = < 47TO' 

l 0, |x| > a, 

in the coordinate system with the origin at the point y, and j(y) is a number 
we can choose as we wish, then c(y) in flMJ) can be calculated: c(y) = r y(y)/30. 
Equation (|23[) implies: 



[V 2 +IC 2 {x)}E e = inM 3 , /C 2 (x) :=if 2 + C(x), (25) 

where C(a;) is defined in (|23p . This equation can be rewritten as 

V x V x E e =)C 2 {x)E e + VV • E e . (26) 

The term VV • E e plays the role of the current iuj/iJ. 

This term can also be interpreted as the term due to a non-local susceptibility 

D e {x) = e(x)E e -iu) I x{ x ,y) E e{y)dy, 



then the Maxwell's equations 

Vx£ e = iio^,H e , V x H e = -iu;e(x)E e - iu x(x,y)E e (y)dy 

Jd 
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imply 

V x V x E e =cj 2 e(x)iiE e (x)+w 2 (i / x{x,y)E e (y)dy 



D 



This equation is of the form J2S|) if e{x) = Jf' , and 

X (x, y) = (w 2 //)" 1 ^^ - y)V y ), 



where S(x — y) is the delta-function. 
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